Improved hourly estimate of aerosol optical thickness over Asian land by fusing geostationary satellites Fengyun-4B and Himawari-9.

Asian over-land aerosols are complexities due to a mixture of anthropogenic air pollutants and natural dust. The accuracy of the aerosol optical thickness (AOT) retrieved from the satellite is crucial to their application in the aerosol data assimilation system. Fusion of AOTs with high spatiotemporal resolution from next-generation geostationary satellites such as Fengyun-4B (FY-4B) and Himawari-9, provides a new high-quality dataset capturing the aerosol spatiotemporal variability for data assimilation. This study develops a complete fusion algorithm to estimate the optimal AOT over-land in Asia from September 2022 to August 2023 at 10 km × 10 km resolution with high efficiency. The data fusion involves four steps: (1) investigating the spatiotemporal variability of FY-4B AOT within the past 1 h and 12 km radius calculation domain; (2) utilizing the aerosol spatiotemporal variability characteristics to estimate FY-4B pure and hourly merged AOTs; (3) performing bias corrections for FY-4B and Himwari-9 hourly merged AOT for different observation times and seasons considering pixel-level errors for each satellite; (4) fusing the bias-corrected FY-4B and Himawari-9 hourly merged AOT based on maximum-likelihood estimation (MLE) method. Compared to the original FY-4B AOT, validation with AERONET observation confirms that the root mean square error (RMSE) of hourly merged FY-4B AOT decreases by around 40.6 % and the correlation coefficient (CORR) increases by about 27.8 %. Compared to FY-4B and Himawari-9 merged AOT, the fused AOT significantly decreases (increases) RMSE (CORR) by around 24.7 % (7.3 %) and 20.2 % (5.6 %). In addition, fused AOT is double the number of single-sensor merged AOT. Fusion aerosol map accurately describes the spatial and temporal variations in Asian regions controlled by air pollution and dust storms. Further studies are required for other landscapes with different satellite combinations to promote the application in the data assimilation system.

Comparison and evaluation of the simulated annual aerosol characteristics over China with two global aerosol models.

In this study, simulations of the annual mean aerosol budget, aerosol optical properties, and surface mass concentration in 2006 in China are performed with two aerosol interactive global atmosphere models, namely, the Nonhydrostatic ICosahedral Atmospheric Model (NICAM) coupled with the Spectral Radiation Transport Model for Aerosol Species (SPRINTARS) and the Beijing Climate Center Atmospheric General Circulation Model (BCC_AGCM) coupled with the Canadian Aerosol Module (CAM) online. The observed and simulated aerosol optical depths (AODs) exhibit similar horizontal distributions with large values over eastern and central China, and sulfate aerosols contribute the main differences between the AODs simulated by NICAM and BCC_AGCM. The simulated sulfate and dust surface concentrations are more consistent with observations compared with the simulated carbonaceous surface concentrations, and both models can reproduce the decreasing tendency of the sulfate surface concentration from urban sites to rural sites. However, the dust emission and deposition levels in China simulated by BCC_AGCM are three times as high as those simulated by NICAM, and the major sink processes of the anthropogenic sulfate, black carbon (BC), and organic carbon (OC) aerosols over China are very different between the two models. The emission and deposition results, which are closely related to the model-assumed aerosol particle size distribution, indicate that the current aerosol size distribution used in the two models should be further improved. The differences in dust emission parameterizations also lead significant discrepancies in aerosol cycles and the dust emission scheme is an important factor determining the magnitudes of global and regional dust emission fluxes.

Reassessment of early 131I inhalation doses by the Fukushima nuclear accident based on atmospheric 137Cs and 131I/137Cs observation data and multi-ensemble of atmospheric transport and deposition models.

The Fukushima Dai-ichi Nuclear Power Plant accidents following the March 11, 2011 Tohoku earthquake, and subsequent tsunami released radioactive materials into the atmosphere and caused significant public health concerns, particularly thyroid cancers in children. However, the lack of measurement data for atmospheric concentrations of 131I has caused persistent and widespread uncertainty. This study estimated the maximum potential thyroid doses of inhaled 131I in the early post-accident phase between March 12 and 23, 2011 by using the hourly measured data of the 137Cs concentrations at 101 suspended particulate matter (SPM) monitoring sites, a new multi-model ensemble (MME) method of simulating 137Cs concentrations using two Atmospheric Transport and Deposition Models (ATDMs), the 131I/137Cs ratio obtained from measurement data analysis, and the internal exposure model. Based on the measurements, the maximum potential thyroid doses were estimated at 3.1-160 mSv at 5 sites in the Fukushima-Hamadori area for 1-year-old children assumed to remain outdoors, whereas they were less than 4.3 mSv at the other sites in the base case of the 131I/137Cs ratio. The spatial distribution of the maximum potential of early inhalation doses was estimated by using the MME and measurements. The inhalation thyroid doses in the evacuation scenarios were compared to the estimates reported by previous studies. The results of the present study were almost congruent with the outcomes of previous investigations except for thyroid doses contributed by highly contaminated plumes on March 12 and 15. The sensitivity analysis for the 131I/137Cs ratio indicated that these plumes carried the potential to significantly increase the thyroid doses of residents.

Effect of global atmospheric aerosol emission change on PM2.5-related health impacts.

Background: Previous research has highlighted the importance of major atmospheric aerosols such as sulfate, through its precursor sulfur dioxide (SO2), black carbon (BC), and organic carbon (OC), and their effect on global climate regimes, specifically on their impact on particulate matter measuring ≤ 2.5 µm (PM2.5). Policy regulations have attempted to address the change in these major active aerosols and their impact on PM2.5, which would presumably have a cascading effect toward the change of health risks. Objective: This study aimed to determine how the change in the global emissions of anthropogenic aerosols affects health, particularly through the change in attributable mortality (AN) and years of life lost (YLL). This study also aimed to explore the importance of using AM/YLL in conveying air pollution health impact message. Methods: The Model for Interdisciplinary Research on Climate was used to estimate the gridded atmospheric PM2.5 by changing the emission of SO2, BC, and OC. Next, the emissions were utilized to estimate the associated cause-specific risks via an integrated exposure-response function, and its consequent health indicators, AM and YLL, per country. Results: OC change yielded the greatest benefit for all country income groups, particularly among low-middle-income countries. Utilizing either AM or YLL did not alter the order of benefits among upper-middle and high-income countries (UMIC/HIC); however, using either health indicator to express the order of benefit varied among low- and low-middle-income countries (LIC/LMIC). Conclusions: Global and country-specific mitigation efforts focusing on OC-related activities would yield substantial health benefits. Substantial aerosol emission reduction would greatly benefit high-emitting countries (i.e. China and India). Although no difference is found in the order of health outcome benefits in UMIC/HIC, caution is warranted in using either AM or YLL for health impact assessment in LIC/LMIC.

Dynamics of atmospheric 131I in radioactive plumes in eastern Japan immediately after the Fukushima accident by analysing published data.

The spatio-temporal distribution of atmospheric radioiodine immediately after the Fukushima Daiichi Nuclear Power Plant (FD1NPP) accident has not yet been clarified due to very limited observed data, compared with atmospheric radiocaesium data. Here, we first revealed that the ratios of 131I (decay-corrected to March 11, 2011) to 137Cs in radioactive plumes were divided into three groups (A, B, and C) by analysing all published data on atmospheric 131I concentrations independently measured immediately after the accident in eastern Japan. Groups A and C were found regardless of whether the measurement sites were located in eastern Fukushima or Kantou areas, while group B was observed only in the eastern Kantou area. The ratios in group A were approximately equal to 10 for the plumes on March 15, March 20, and on the morning of March 21, and those in group B were approximately 75 on March 16. Their possible sources were Unit 2 and/or Unit 3. In contrast, the ratios in group C were approximately equal to 360, much higher than those of groups A and B, and were observed from the afternoon of March 21 to March 25. These high 131I concentrations could be released after water supply to FD1NPP.

Aerosol effects on cloud water amounts were successfully simulated by a global cloud-system resolving model.

Aerosols affect climate by modifying cloud properties through their role as cloud condensation nuclei or ice nuclei, called aerosol-cloud interactions. In most global climate models (GCMs), the aerosol-cloud interactions are represented by empirical parameterisations, in which the mass of cloud liquid water (LWP) is assumed to increase monotonically with increasing aerosol loading. Recent satellite observations, however, have yielded contradictory results: LWP can decrease with increasing aerosol loading. This difference implies that GCMs overestimate the aerosol effect, but the reasons for the difference are not obvious. Here, we reproduce satellite-observed LWP responses using a global simulation with explicit representations of cloud microphysics, instead of the parameterisations. Our analyses reveal that the decrease in LWP originates from the response of evaporation and condensation processes to aerosol perturbations, which are not represented in GCMs. The explicit representation of cloud microphysics in global scale modelling reduces the uncertainty of climate prediction.

Estimation of chlorophyll concentration in waters near Hokkaido using the linear combination method.

An inversion algorithm is implemented to retrieve the surface chlorophyll a (Chl) concentration using satellite observation data from the MODIS instrument. The algorithm employs a simple and flexible index (LCI) to combine with the Chl without explicit correction for aerosol scattering. To investigate the sensitivity of LCI to Chl and other influence factors, an oceanic radiative transfer model coupled with a comprehensive bio-optical module is developed. It is studied that the LCI is significantly linear to Chl and not sensitive to other influence factors, except in very low oceanic salinity or scattering angle conditions, where over a 12% relative difference of derived Chl exists. Inversion results show the retrieved Chl are highly consistent with the MODIS operational data products in waters near Hokkaido, with the correlation coefficient, root mean square deviation, and average percentage difference of 0.9702, 0.3756 mg m-3, and 13.89%, respectively. Investigation of the validity of this algorithm with a variety of atmospheric conditions indicates that the residual influence of atmosphere on the LCI index, after Rayleigh scattering correction, is generally within ± 0.001, allowing the retrieval error of Chl at less than 25% in most cases. A good comparison between retrieval and in situ measurements is also identified and implies that the retrieval accuracy via the LCI method depends on the linear combination coefficients used and the bio-optical module selected, while effects of polarization can be ignored.

Pseudoaneurysm formation due to rupture of intracranial aneurysms: Case series and literature review.

Background Intracranial pseudoaneurysm formation due to a ruptured non-traumatic aneurysm is extremely rare. We describe the radiological findings and management of pseudoaneurysms due to ruptured cerebral aneurysms in our case series and previously reported cases. Patients and methods Four additional and 20 reported patients presenting with subarachnoid hemorrhage (SAH) are included. Radiological findings and clinical features of these patients were reviewed. Results In our series, three-dimensional computed tomographic angiography (3D-CTA) and/or angiography showed an irregular- or snowman-shaped cavity extending from the parent artery. The radiological examination additionally revealed delayed filling and retention of contrast medium. These findings were the same as previously reported cases. One patient underwent direct clipping of the true aneurysm. For the other three patients with aneurysms at the basilar and anterior communicating arteries, the true portion of the aneurysm was embolized with platinum coils. During the procedures, care was taken not to insert the coils into the distal pseudoaneurysm portion to prevent rupture. The review of 24 cases revealed that the location of the aneurysms was most frequent in the anterior communicating artery (41.7%), and 86.7% of patients were in a severe stage of SAH (>Grade 3 in WFNS or Hunt & Kosnik grading) implying abundant SAH. Conclusions Pseudoaneurysm formation in SAH after non-traumatic aneurysm rupture is rare. However, in cases with an irregular-shaped aneurysm cavity, pseudoaneurysm formation should be taken into consideration.

Reducing Multi-sensor Monthly Mean Aerosol Optical Depth Uncertainty Part II: Optimal Locations for Potential Ground Observation Deployments.

Surface remote sensing of aerosol properties provides "ground truth" for satellite and model validation, and is an important component of aerosol observation system. Due to the different characteristics of background aerosol variability, information obtained at different locations usually have different spatial representativeness, implying that the location should be carefully chosen so that its measurement could be extended to a greater area. In this study, we present an objective observation array design technique that automatically determines the optimal locations with the highest spatial representativeness based on the Ensemble Kalman Filter (EnKF) theory. The ensemble is constructed using aerosol optical depth (AOD) products from five satellite sensors. The optimal locations are solved sequentially by minimizing the total analysis error variance, which means that observations at these locations will reduce the background error variance to the largest extent. The location determined by the algorithm is further verified to have larger spatial representativeness than some other arbitrary location. In addition to the existing active AERONET sites, the 40 selected optimal locations are mostly concentrated on regions with both high AOD inhomogeneity and its spatial representativeness, namely the Sahel, South Africa, East Asia and North Pacific Islands. These places should be the focuses of establishing future AERONET sites in order to further reduce the uncertainty in the monthly mean AOD. Observations at these locations contribute to approximately 50% of the total background uncertainty reduction.

Inverse modeling of the 137Cs source term of the Fukushima Dai-ichi Nuclear Power Plant accident constrained by a deposition map monitored by aircraft.

The amount of 137Cs released by the Fukushima Dai-ichi Nuclear Power Plant accident of 11 March 2011 was inversely estimated by integrating an atmospheric dispersion model, an a priori source term, and map of deposition recorded by aircraft. An a posteriori source term refined finer (hourly) variations comparing with the a priori term, and estimated 137Cs released 11 March to 2 April to be 8.12 PBq. Although time series of the a posteriori source term was generally similar to those of the a priori source term, notable modifications were found in the periods when the a posteriori source term was well-constrained by the observations. Spatial pattern of 137Cs deposition with the a posteriori source term showed better agreement with the 137Cs deposition monitored by aircraft. The a posteriori source term increased 137Cs deposition in the Naka-dori region (the central part of Fukushima Prefecture) by 32.9%, and considerably improved the underestimated a priori 137Cs deposition. Observed values of deposition measured at 16 stations and surface atmospheric concentrations collected on a filter tape of suspended particulate matter were used for validation of the a posteriori results. A great improvement was found in surface atmospheric concentration on 15 March; the a posteriori source term reduced root mean square error, normalized mean error, and normalized mean bias by 13.4, 22.3, and 92.0% for the hourly values, respectively. However, limited improvements were observed in some periods and areas due to the difficulty in simulating accurate wind fields and the lack of the observational constraints.

Unrealistically pristine air in the Arctic produced by current global scale models.

Black carbon aerosol (BCA) in the Arctic has profound impacts on the global climate system through radiation processes. Despite its enormous impacts, current global scale models, powerful tools for estimating overall impact, tend to underestimate the levels of BCA in the Arctic over several seasons. Using a global aerosol transport simulation with a horizontal grid resolution of 3.5 km, we determined that a higher resolution significantly reduced the underestimation of BCA levels in the Arctic, mainly due to an enhancement of the representation of low-pressure and frontal systems. The BCA mass loading in the Arctic simulated with 3.5-km grid resolution was 4.2-times larger than that simulated with coarse (56-km) grid resolution. Our results also indicated that grid convergence had not occurred on both the contrast between the cloud/cloud free areas and the poleward BCA mass flux, despite the use of the 3.5-km grid resolution. These results suggest that a global aerosol transport simulation using kilometre-order or finer grid resolution is required for more accurate estimation of the distribution of pollutants in the Arctic.

Reducing multisensor satellite monthly mean aerosol optical depth uncertainty: 1. Objective assessment of current AERONET locations.

Various space-based sensors have been designed and corresponding algorithms developed to retrieve aerosol optical depth (AOD), the very basic aerosol optical property, yet considerable disagreement still exists across these different satellite data sets. Surface-based observations aim to provide ground truth for validating satellite data; hence, their deployment locations should preferably contain as much spatial information as possible, i.e., high spatial representativeness. Using a novel Ensemble Kalman Filter (EnKF)-based approach, we objectively evaluate the spatial representativeness of current Aerosol Robotic Network (AERONET) sites. Multisensor monthly mean AOD data sets from Moderate Resolution Imaging Spectroradiometer, Multiangle Imaging Spectroradiometer, Sea-viewing Wide Field-of-view Sensor, Ozone Monitoring Instrument, and Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar are combined into a 605-member ensemble, and AERONET data are considered as the observations to be assimilated into this ensemble using the EnKF. The assessment is made by comparing the analysis error variance (that has been constrained by ground-based measurements), with the background error variance (based on satellite data alone). Results show that the total uncertainty is reduced by ~27% on average and could reach above 50% over certain places. The uncertainty reduction pattern also has distinct seasonal patterns, corresponding to the spatial distribution of seasonally varying aerosol types, such as dust in the spring for Northern Hemisphere and biomass burning in the fall for Southern Hemisphere. Dust and biomass burning sites have the highest spatial representativeness, rural and oceanic sites can also represent moderate spatial information, whereas the representativeness of urban sites is relatively localized. A spatial score ranging from 1 to 3 is assigned to each AERONET site based on the uncertainty reduction, indicating its representativeness level.

First retrieval of hourly atmospheric radionuclides just after the Fukushima accident by analyzing filter-tapes of operational air pollution monitoring stations.

No observed data have been found in the Fukushima Prefecture (FP) for the time-series of atmospheric radionuclides concentrations just after the Fukushima Daiichi Nuclear Power Plant (FD1NPP) accident. Accordingly, current estimates of internal radiation doses from inhalation, and atmospheric radionuclide concentrations by atmospheric transport models are highly uncertain. Here, we present a new method for retrieving the hourly atmospheric (137)Cs concentrations by measuring the radioactivity of suspended particulate matter (SPM) collected on filter tapes in SPM monitors which were operated even after the accident. This new dataset focused on the period of March 12-23, 2011 just after the accident, when massive radioactive materials were released from the FD1NPP to the atmosphere. Overall, 40 sites of the more than 400 sites in the air quality monitoring stations in eastern Japan were studied. For the first time, we show the spatio-temporal variation of atmospheric (137)Cs concentrations in the FP and the Tokyo Metropolitan Area (TMA) located more than 170 km southwest of the FD1NPP. The comprehensive dataset revealed how the polluted air masses were transported to the FP and TMA, and can be used to re-evaluate internal exposure, time-series radionuclides release rates, and atmospheric transport models.

Improvement of aerosol optical properties modeling over Eastern Asia with MODIS AOD assimilation in a global non-hydrostatic icosahedral aerosol transport model.

A new global aerosol assimilation system adopting a more complex icosahedral grid configuration is developed. Sensitivity tests for the assimilation system are performed utilizing satellite retrieved aerosol optical depth (AOD) from the Moderate Resolution Imaging Spectroradiometer (MODIS), and the results over Eastern Asia are analyzed. The assimilated results are validated through independent Aerosol Robotic Network (AERONET) observations. Our results reveal that the ensemble and local patch sizes have little effect on the assimilation performance, whereas the ensemble perturbation method has the largest effect. Assimilation leads to significantly positive effect on the simulated AOD field, improving agreement with all of the 12 AERONET sites over the Eastern Asia based on both the correlation coefficient and the root mean square difference (assimilation efficiency). Meanwhile, better agreement of the Ångström Exponent (AE) field is achieved for 8 of the 12 sites due to the assimilation of AOD only.

Development of a light scattering solver applicable to particles of arbitrary shape on the basis of the surface-integral equations method of Müller type. I. Methodology, accuracy of calculation, and electromagnetic current on the particle surface.

We develop a numerical algorithm for calculating the light-scattering properties of small particles of arbitrary shape on the basis of a method involving surface integral equations. The calculation error was estimated by performing a comparison between the proposed method and the exact Mie method with regard to the extinction efficiency factor, and the results show that the error is less than 1% when four or more nodes per wavelength are set on the surface of a spherical particle. The accuracy fluctuates in accordance with the distribution of nodal points on the particle surface with respect to the direction of propagation of the incident light. From our examinations, it is shown that the polar incidence alignment yields higher accuracy than equator incidence when a "latitude-longitude" type of mesh generation is adopted. The electric currents on the particle surface and the phase functions of all scattering directions are shown for particles shaped as spheres or hexagonal columns. It is shown that the phase function for a hexagonal column has four or eight cold spots. The phase function of a randomly oriented hexagonal column shows halolike peaks with size parameters of up to 20. This method can be applied to particles with a size parameter of up to about 20 without using the symmetry characteristic of the particle.

Application of the SKYRAD Improved Langley plot method for the in situ calibration of CIMEL Sun-sky photometers.

The in situ procedure for determining the solar calibration constants, originally developed for the PREDE Sun-sky radiometers and based on a modified version of the Langley plot, was applied to a CIMEL instrument located in Valencia, Spain, not integrated into AERONET. Taking into account the different mechanical and electronic characteristics of the two radiometers, the method was adapted to the characteristics of the CIMEL instrument. The iterative procedure for the determination of the solar calibration constants was applied to a 3-year data set. The results were compared with the two sets of experimental calibration constants determined during this period using the standard Langley plot method. The agreement was found to be consistent with the experimental errors, and the method can definitely also be used to determine the solar calibration constant for the CIMEL instrument, improving its calibration. The method can be used provided the radiometer is previously calibrated for diffuse radiance using a standard lamp.

Determination of the solar calibration constant for a sun-sky radiometer: proposal of an in-situ procedure.

A new procedure is presented for determining in situ the solar calibration constant, i.e., the Sun-sky radiometer counts for a direct normal solar flux extrapolated to the top of the atmosphere. The method makes use of a modified version of the Langley plot based on the use of an inversion code of column-integrated aerosol size distribution, and it is ordinarily applied to calibrate Prede Sun-sky radiometers. To analyze how such an in situ method can work accurately, the technique has been applied to a five-month dataset obtained from measurements taken in Rome, Italy, by a Prede Sun-sky radiometer from 22 April to 5 November 2001. The precision of the in situ method has been estimated to within 1-2.5%, depending on the wavelength.

Efficient use of an improved radiative transfer code to simulate near-global distributions of satellite-measured radiances.

Two new extension modules that give the water-leaving radiance from the ocean and the snow bidirectional reflectance distribution function were implemented in the latest radiative transfer code. In addition, to simulate the near-global distributions of satellite-measured radiances by using the improved radiative transfer code, we tested and applied the look-up table method together with the process-separation technique of the radiative transfer calculation. The computing time was reduced from 1 year to 20 s to simulate one channel, one scene of the Global Imager image by use of an Alpha 21164A-2 (600-MHz) machine. The error analyses showed that the radiances were simulated with less than 1% error for the nonabsorbing visible channels and approximately 2% error for absorbing channels by use of this method.